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The influence of deviations in dwell times and source positions for ^'^Ir HDR-RALS was investigated. The 
potential dose errors for various kinds of brachytherapy procedures were evaluated. The deviations of dwell 
timeArofa^'^^IrHDR source for the various dwell times were measured with a well-type ionization chamber. 
The deviations of source position AP were measured with two methods. One is to measure actual source pos- 
ition using a check ruler device. The other is to analyze peak distances from radiographic film irradiated with 
20 mm gap between the dwell positions. The composite dose errors were calculated using Gaussian distribu- 
tion with Arand AP as lo of the measurements. Dose errors depend on dwell time and distance from the point 
of interest to the dwell position. To evaluate the dose error in clinical practice, dwell times and point of interest 
distances were obtained from actual treatment plans involving cylinder, tandem-ovoid, tandem-ovoid with 
interstitial needles, multiple interstitial needles, and surface-mold applicators. The Ar and AP were 32 ms 
(maximum for various dwell times) and 0.12 mm (ruler), 0.11 mm (radiographic film). The multiple intersti- 
tial needles represent the highest dose error of 2%, while the others represent less than approximately 1%. 
Potential dose error due to dwell time and source position deviation can depend on kinds of brachytherapy 
techniques. In all cases, the multiple interstitial needles is most susceptible. 
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INTRODUCTION 

The radioactive source and seed used in brachytherapy can emit 
photons or electrons in all directions. DosimetricaUy, this means 
that the dose gradient around a brachytherapy source is very 
steep, and the dose decreases exponentially with the thickness 
of the material or tissue. Therefore, by implanting the source 
near or into the tumor region, high conformity of prescribed 
dose to the target and reduction of unwanted dose to healthy 
organs surrounding the tumor region can be obtained. 

The use of CT or MR images in treatment planning for 
brachytherapy has been widely employed as image-guided 
brachytherapy (IGBT). Over the past few years, a number of 
studies have been conducted on various IGBT techniques [1^]. 
Dose distributions can be customized to fit the target for each 
patient by referring to actual anatomy from CT or MR 
images instead of film-based planning. Furthermore, recently, 
sophisticated dwell time optimization techniques, such as 
graphical optimization and inverse planning anatomy-based 



dose optimization, has been implemented in commercial treat- 
ment planning systems (TPSs) [5-9]. 

Treatment parameters (including dwell times and dwell 
positions of the radioactive source calculated by TPSs) are 
transferred to the treatment machine, such as a remote after- 
loading system (RALS), manufactured e.g. by Nucletron, an 
Elekta company (MicroSelectron®, Stockholm, Sweden), 
Varian Medical Systems Inc. (VariSource™, Palo Alto, CA, 
USA), or Eckert & Ziegler Bebig SA (MultiSource®, Seneffe, 
Belgium). The RALS should provide high precision control of 
the movement of the radioactive source as planned in terms of 
the dwell times and dwell positions in the applicators in 
which the source moves. Mechanical accuracies of move- 
ment of the radioactive source have been referred to in some 
quality assurance/quality control (QA/QC) guidelines for 
high-dose-rate (HDR) RALS [10-13]. For example, AAPM 
Task Group 40 and Task Group 56 refer to a tolerance of 
source positioning of ±1 mm [10, 11]. ESTRO Booklet No. 8 
recommend an action level of source positioning accuracy 
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of ±2 mm [12]. The tolerance or action level can be applic- 
able for the scheduled QA/QC treatment machine. However, 
it is not obvious how these figures relate to clinical practice, 
and they lack practical value in terms of clinical influence. It 
has been further suggested that these values can be applied 
to all of the irradiation techniques in brachytherapy, such as 
surface-mold and interstitial brachytherapy, etc. Actually, 
dose error due to deviation of the source position depends on 
the distance from the source according to the inverse square 
law. A large dose error can be observed near the source, but 
the dose error far from the source is less. 

Similarly, the impact of dose change caused by dwell time 
error can be expected to be linearly dependent on total dwell 
time. For instance, a tiny dwell time error can have a large 
impact in the case of low dwell time. In particular, the graph- 
ical optimization and inverse-planning anatomy-based dose 
optimization of the commercially available TPSs performs 
calculations of the dwell times without putting a lower limit 
on dwell time. Therefore, it is important to investigate dose 
error for such brachytherapy techniques, which freely use a 
low dwell time. 

Dose error from mechanical uncertainties of movement of 
the source can depend on the kind of irradiation technique used, 
and its clinical influence is not fully understood. Our purpose in 
this study was to evaluate the dose error for clinical treatment 
plans involving a range of irradiation apphcators/catheters, such 
as cylinder, tandem/ovoid, tandem/ovoid with a small number 
of interstitial needle applicators ('Combination-brachytherapy'), 
multiple interstitial needle applicators, and surface-mold appli- 
cators for scalp tumors. 

In this study, deviations in dwell time and source position 
were measured for a number of techniques. Finally, the influ- 
ence of the dose error due to both deviations was investigated 
for actual treatment plans involving the various applicators 
or catheters listed above. 



Dwell time can be obtained from the amount of net 
electric charge (M - Mef) and the current Aq of the WIC 
during the irradiation, as shown in the following equation. 

M - Mef /,^ 
Tm = ■ (1) 

Mef means amount of electric charge in the end effect of a 
source, and it can be regarded as a constant independent of 
dwell time. Therefore, a standard deviation of dwell time T,„ 
can be equal to that of M/Aq. The dwell time deviation AT 
was defined as a maximum value for a standard deviation of 
10 measurements of MIAq for nominal dwell time T of 0.1, 
0.2, 0.5, 1.0,5.0, 15.0 and 30.0 s. 

The deviations in the source positions were measured with 
two methods. One method involved measuring the position 
of the source with a source position check ruler device spe- 
cifically developed for MicroSelectron® and provided by 
Nucletron, an Elekta company. The MicroSelectron® was set 
up to place the source in a certain position via the ruler con- 
nected by a straight transfer tube. By reading the actual 
source position, the standard deviation of the position devi- 
ation was obtained. A total of 11 measurements were per- 
formed to obtain a standard deviation. The other method 
involved irradiation of radiographic film (EDR2, Carestream 
Health Inc., Rochester, NY, USA) with applicators tightly 
attached to it. Irradiation was performed with 2-cm gaps 
between the dwell positions, and the number of gaps was 
198. The exposed films were digitized with an EPSON 
ES-8500 flatbed scanner under resolution conditions of 
600 dpi (equivalent to 0.042 mm/pixel) and a 16-bit gray 
scale. The digitized images were analyzed using film QA soft- 
ware (DD-system Ver. 10.21, R-TECH Company, Tokyo, 
Japan) as shown in Fig. 1. The source position deviation AP 
was defined as the maximum of a standard deviation for two 
tests (ruler and radiographic film). 



MATERIALS AND METHODS 

Measurements of deviation of dwell time and source 
position 

For ''^Ir HDR-RALS (MicroSelectron® V2, Nucletron, an 
Elekta company), the deviation in dwell time AT and source 
position AP were measured for a range of techniques [12-16]. 
For two measurements, the air kerma stiiength was -30.0 mGy 
m^ h~'. The Nucletron RALS has a nominal time resolution of 
0.1 s. As shown in Eq. (1), dwell time, T„, was obtained from 
the amount of electric charge measured by a well-type ioniza- 
tion chamber (WIC) and an electrometer (HDR-1000 plus and 
MAX-4000, Standard Imaging Inc., Middleton, WI, USA). 
The correction factor kxp was applied to convert the cavity air 
mass at the reference conditions (temperature Tq 22°C and 
pressurePo 101.3 kPa) [13]. 



Evaluation of dose error in single catheter method 

As shown in Fig. 2, we calculated the dose errors at the 
points of interest due to deviations in both dwell time AT and 
source position AP, using the methods as described in the 
previous section. Two deviations were used as a standard de- 
viation (according to Gaussian distribution) in the develop- 
ment environment of Microsoft Visual Studio 2010 instead 
of the commercial TPS. Dose calculation with a grid size of 
1 .0 mm was performed according to AAPM TG-43 protocol 
[17-19], as is typically used for dose calculation in brachy- 
therapy TPS. Points of interest were placed at a distance L 
from each dwell position in the perpendicular direction of a 
catheter. A total of 10 dwell positions with a 2.5-mm dis- 
tance from each source was investigated. Dose deviations 
were defined as the mean value for a standard deviation of 
the dose at the points of interest in 10 trials. 
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Fig. 2. Dose errors at points of interest for a single catheter 
calculated according to the AAPM TG-43 protocol [17-19] due to 
the two deviations: dwell time deviation Ar and source position 
deviation AP. 



Evaluation of dose calculation error in various 
clinical treatment plans 

Dose error at the points of interest is affected by two factors: 
dwell time and distance of the point of interest from the source. 
In order to evaluate the influence of the two deviations, zirand 
AP, for a clinical treatment plan, we measured the dwell times 
and distance to points of interest in the actual treatment plan 
for a range of applicator types, such as cylinder, tandem/ovoid, 
combination-brachytherapy, multiple interstitial needle appli- 
cators, and surface-mold applicators for scalp tumors, in TPS 
Oncentra Brachy Ver. 4.1.0.132 (Nucletron, Veenendaal, 
Netherland). Table 1 provides information about the treatment 
plans for the various applicator types investigated here. 

As time passes, dwell times for the same treatment plan 
must be increased in order to compensate for the decay of the 
radioactivity source (half-Ufe of '''^Ir = 73.83 d). In other words, 
air kerma strength is higher, e.g., just after the replacement of an 
^^^Ir source. A maximum value Sk of 48.56 mGy m^ h~^ was 
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Applicators or catheters 


# treatnient 
plans 


Prescribed dose 
(Gy) 


Method of optimizing dose distribution 


Cylinder 


16 


6 


Dose point optimization to 5 mm under vaginal surface 


Tandem/ovoid 


17 


6 


Manchester system 


Combination-brachytherapy 


9 


6 


Inverse planning based on dose-volume histogram 
and graphical optimization 


Multiple interstitial needle 
applicators 


8 


6 




Surface-mold for scalp tumor 


5 


2-2.5 


Dose point and graphical optimization 



obtained from the 43 measurements taken between February 
2002 and December 2012 in our institute (average ± ISD, 
42.21 ± 3.06 mGy h~^). Dwell times for all of the treatment 
plans hsted in Table 1 were recalculated using the maximum 
Sfc of 48.56 mGy m^ h~\ because a large dose error can be 
observed in the case of low dwell time. 

In this study, the distance to the point of interest was 
defined by the equation below. We used the geometrical 
averaged distance of the point of interest from catheter, r^. 




This equation can be derived from ^ ~ ^ioo%- 

V\oo% means the irradiated volume redeiving the prescri- 
bed dose or more, and it can be obtained from the dose- 
volume histogram in the TPS. The /, is assumed to be the 
length of treatment region of the ith catheter. The gap 
between the dwell positions was set to a 2.5-mm interval. 
Therefore, the length of treatment region /, could be 
expressed as a multiplication of the 2.5 mm interval by the 
number of dwell positions. 

RESULTS 

Measurements of deviation of dwell time 
and source position 

Dwell time deviation was measured with a WIC. The stand- 
ard deviations of dwell time for 10 measurements are shown 
in Table 2. From these results, the maximum standard devi- 
ation was found to be 32 ms for a dwell time of 30.0 s. 

Deviation of the source position in the Nucletron RALS 
was measured using the two methods. The standard deviation 
of the source position was 0.12 mm and 0.11 mm in mea- 
surements using either the source position check ruler device 
(«=11) or radiographic film (m=198), respectively. The 
dwell time and source position deviation in this study were 
regarded as a standard deviation according to Gaussian distri- 
bution, AT= 32 ms and AP = 0. 12 mm (worst case scenario). 



Table 2. Dwell time deviation for Nucletron RALS 
measured with a WIC 

Dwell time (s) 0.1 0.2 0.5 1.0 5.0 15.0 30.0 
Deviation (ms) 18 17 23 23 29 23 32 



Dose error in a single catheter by calculation 

Figure 3a depicts the calculated dose errors as a function of 
dwell time with points of interest placed at a fixed distance of 
5 mm with the measured two deviations, using AT= 32 ms 
and AP = 0.12 rrtm as the standard deviation according to the 
Gaussian distribution. Dose errors were calculated with three 
conditions: with either AT, AP, or both of them to reveal the 
contribution of each deviation to the total dose error. Using 
only AP, the dose error was constant independent of the 
dwell time. Using only AT, the dose error decreased with an 
increased dwell time. For both two deviations, the dose error 
also decreased with an increase in dwell time, and was a com- 
posite of the graphs when using only AP and only AT. Figure 
3b depicts the calculated dose errors as a function of distance 
L to the point of interest with a fixed dwell time of 2 s and the 
two measured deviations, AT= 32 ms and AP = 0.12 mm gen- 
erating the standard deviation according to Gaussian distribu- 
tion. Dose errors were also calculated according to these 
assumptions in Fig. 3a. With only zlP, the dose error decreased 
with an increase in distance L to the point of interest. The dose 
error approached zero, because zlP is less influential for a 
point distant from the source. Using only AT, the dose error 
remained constant, because dose was linearly proportional to 
the dwell time. With both deviations, the dose error also 
decreased with increase in distance to the point of interest, and 
was a composite of the graphs with only AP and only AP. 

As shown in Fig. 3, the dose error depended on two factors: 
dwell time and distance to the point of interest. Figure 4 
depicts a 2D-plot of the dose errors with the two measured 
deviations, Ar=32ms and AP = 0.12mm. As shown in the 
figure, a high dose error could be observed near the source for 
a small dwell time. 
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Fig. 3. Dose error for Ar=32ms and A/' = 0.12mm: (a) dose 
error at points of interest placed at a fixed distance of 5 mm as a 
function of dwell time by calculation, (b) Dose error for a fixed 
dwell time of 2 s as a function of distance to points of interest by 
calculation. 



Dwell time and distance to the point of interest 
for clinical treatment plans 

Figure 5 depicts the box-and-whisker plot of dwell times per 
dwell position in the actual clinical treatment plans as 
described in Table 1. The median dwell time for each dwell 
position was 27.18, 19.10, 5.74, 0.73 and 0.95 s for the cylin- 
der, tandem/ovoid, combination-brachytherapy, multiple inter- 
stitial needle appUcators, and surface-mold applicators for 
scalp tumors, respectively. The lower quartile of dwell times 
was 20.96, 13.26, 1.97, 0.34 and 0.69 s, respectively. For in- 
stance. Fig. 6 depicts a distribution of dwell times for the mul- 
tiple interstitial needle appUcators (n = 8 treatment plans). 

Figure 7 depicts the box-and-whisker plot of the geometrical 
averaged distance to the point of interest derived from Eq. 
(2) for the clinical treatment plans (see Table 1). The median 
value of rg was 30.73, 29.06, 17.01, 8.08 and 20.10 mm in the 
cylinder, tandem/ovoid, combination-brachytherapy, multiple 
interstitial needle applicators, and surface-mold applicators for 
scalp tumors, respectively. The lower quartile of was 30.14, 
28.51, 16.37, 7.44 and 15.70 mm, respectively. Furthermore, 




1 2 3 4 5 10 15 20 25 30 
Distance to point of interest L (mm) 



Fig. 4. 2D-plot of dose error for dwell times and distance to point 
of interest L for Ar = 32 ms and hP = 0.12 mm. 
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Fig. 5. The box-and-whisker plot of the dwell times per source 
position for the clinical treatment plans (see Table 1). 



minimum rg was 27.03, 27.22, 14.72, 6.47, and 15.29 mm, 
respectively. 

Potential dose error for clinical treatment plans 

We adopted the results of the dose error for a single catheter 
by calculations from clinical treatment plans. We introduced 
the median value of Figs 5 and 7 as a representative expres- 
sion of the treatment plans. This condition assumed that the 
dwell time and the distance to the point of interest for all acti- 
vations can be considered equivalent to the median value. 
Table 3 lists the potential dose errors that were obtained from 
the 2D-plot of the dose error indicated by the two median 
values. 
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Fig. 6. Distribution of dwell times for the multiple interstitial 
needle applicators (n = 8 treatment plans). Inverse planning and 
graphical optimization are used to optimize dwell time in the 
Treatment Planning System. 
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Fig. 7. The box-and-whisker plot of the geometrical averaged 
distance to the point of interest derived from Eq. (2) for clinical 
treatment plans (see Table 1). 



Table 3. Potential dose errors for clinical treatment plans 
(see Table 1) with Ar=32ms and AP = 0. 12 mm from 
measurements 



Applicators or catheters 



{AT, AP) = (32 ms, 0.12 mm) 
Potential dose error 



Cylinder 


-0.1% 


Tandem/ovoid 


-0.2% 


Combination-brachytherapy 


-0.3% 


Multiple interstitial needle 


-2.0% 


applicators 




Surface-mold for scalp tumor 


-1.3% 


DISCUSSION 



In brachytherapy, there are many uncertainties in the pre- 
scribed dose for the overall treatment process such as the 



uncertainty in calculating the absolute dose [17-19] and 
the effect of inhomogeneity [20-24], the mechanical ac- 
curacy of the source movement, and the intra- and inter- 
fractional displacements of the applicators or catheters, 
which cause substantial dose deviation from the planned 
dose. Of these uncertainties, we focus upon the dose errors 
caused only by the mechanical uncertainties of movement of 
the source in '''"Ir HDR-RALS in terms of dwell time and 
source position. 

The dwell time deviation was obtained from the amount 
of electric charge from the WIC during the irradiation [see 
Eq. (I)]. This measurement is considered to be influenced 
by the deviation, both in dwell time and source position. 
However, we could neglect the influence of deviation in the 
source position, for the following reason. For this measure- 
ment, the ^'^Ir souire was located at a predetermined position 
in the WIC where the amount of electric charge represented 
the maximum. The behavior of the amount of electric charge 
in this region can be comparatively flat for a variable source 
position. For example, a displacement of 0. 1 mm in the source 
position causes a change in the amount of electric charge by 
0.001%. This corresponds to a dwell time of 1 ms. Therefore, 
the amount of electric charge is weakly dependent on the devi- 
ation of the source position. 

The Nucletron RALS has high precision control of the 
dwell time of a source independent of dwell time (Table 2). 
We determined source position deviation via two methods, 
and the results obtained from each in this study were similcU". 
Mechanical accuracy of the source movement in the short 
term was also evaluated. Actually, deterioration of the motor 
controlling source movement can happen abruptly or can 
gradually worsen in the long term. The long-term mechanic- 
al accuracy needs to be further investigated. 

From Fig. 3, it is clear how dose eiTor at the point of inter- 
est for a single catheter can change with variable dwell time 
and distance to the point of interest from a source. In the case 
of low dwell time and placement of the point of interest near 
a source, dose error can be expected to be high. Therefore, it 
is important that the effect of both deviations is taken into 
account when evaluating potential dose error for different the 
brachytherapy techniques. Of all the brachytherapy tech- 
niques investigated, the smallest and largest dwell times 
were obtained with the multiple interstitial needle and cylin- 
der applicators, respectively (see Fig. 5). Additionally, the 
brachytherapy techniques having the smallest and largest dis- 
tance to the point of interest from the source were also the 
multiple interstitial needle and cylinder applicators, respect- 
ively (see Fig. 7). 

As shown in Table 3, all of the catheters or applicators 
represented fell within a 2% dose error. This is because the 
Nucletron RALS has high precision control of the move- 
ment of the radioactive '''"Ir source, Ar=32ms and 
zip = 0.12 mm. The finding is that the potential dose error 
depends on the kind of brachytherapy technique used. Of all 
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the techniques investigated, the multiple interstitial needle 
applicators could introduce the most error, -2.0% potential 
dose error. A procedure for evaluation of the potential dose 
error, applicable for all brachy therapy techniques, was pro- 
posed in this study. The authors recommended that the user 
establish the 2D-plot for dwell time and distance to the 
point of interest with the measured accuracy of movement 
of the source, and need to evaluate the potential dose error 
caused by mechanical uncertainties of movement of the 
source for the brachytherapy techniques used in a particular 
institute. 

CONCLUSIONS 

Mechanical uncertainties in the '^^Ir HDR-RALS (Nucletron, 
an Elekta company), a dwell time deviation AT of 32 ms and a 
source position deviation AP of 0.12 mm were obtained from 
measurements. The finding was that the Nucletron HDR-RALS 
has high precision control of the movement of the radioactive 
^'^Ir source in terms of dwell time and source position. The 
potential dose error caused by the two deviations depends on 
the kind of brachytherapy technique used. Of the techni- 
ques studied, multiple interstitial needle applicators have 
the highest susceptibility to dwell time deviation and source 
position deviation. 
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